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ANAIiYTICAL 2U® NUMERICAL INVESTIGATION OP STRUCTURAL 
RESPONSE OP COMPLIANT WALL MATERIALS 

By 

R. Balasubramanicin^ 


SUMMARY 


Theoretical analysis of an electrostatically driven wall system for a 
Gomp3j.ant wall drag reduction program is reported. The electrostatic wall 
system is capable of producing deflections of many orders greater than the 
thicknesses and at small wavelengths. An intermediate large response theory 
is used for structural analysis. The theoretical predictions are compared to 
bench test results, and good agreement between the two is obtained. The effects 
of aerodynamic loads and perturbation electric fields on the theoretical 
solutions are considered. It is shown that for very small wavelengths 
(X « 2 mm) the aerodynamic effects can be estimated using potential theory 
without loss of accuracy, and the perturbation electric fields do not affect 
solutions as long as the deflections are less than one percent of the wave- 
length. Resonance effects for this type of structure are shown to be fairly 
small. 


1. INTRODUCTION 


Details of the compliant wall drag reduction program at Langley have 
been discussed in a supplementary report (ref. 1) under the present contract 
NSG 1236. It was pointed out in that report that passive walls with short 
wavelengths and large amplitudes are extremely difficult to design. Extension 
of grant NSG 1236 was given in order to design controlled active wall experi- 
ments. The amplitude of surface motion desired was given to be in the range 
of 5 X 10“^ m to 2 X 10“*^ m, the wavelength of the s\n:face motion in the range 
2,5 X 10“^ m to 6 X 10“^ m, and the frequency range to be between 300 Hz and 
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2 kHz. After careful evaluation of existing techniques for active wall 
experiments (refs. 2 to 4) it was decided to develop an active wall system using 
electrostatic forces as loading to the structure. The choice of materials 
for this systKti was narrowed to extremely thin elastomers, and the amplitude 
constraint given above dictated large values of the ratio of amplitude thick- 
ness. A nonlinear structural response analysis was conducted to determine 
accurate surface motion predictions, in sections 2 to 4 the electrostatic 
wall system is discussed in detail, and in section 5 comparisons with e:q)eri- 
mental measurements are reported. 

2. ANALYSIS OF THE ELECTROSTATIC WALL PROBLEM 

The recent summary paper (ref. 5) suggests that low-speed air experiments 
for compliant wall drag reduction be conducted with controlled or active wall sur- 
face motion to assess the nature of possible turbulent boundary layer modifications 
due to the wall motion. Previous active wall experiments reported in the litera- 
ture (refs. 2 to 4) used mechanical drivers; these drive systems are inadequate for 
producing high-frequency, short-wavelength motion contemplated for active wall experi- 
ments at Leuigley. The present work describes an electrostatic wall designed to 
operate in a frequency range of 200 Hz to 10 kHz with two-dimensional standing waves 
of wavelength 2 x 10“^ to 10“^ m. The structural surface is basically a thin elec- 
trically conducting elastomer membrane with a series of transverse electrodes etched 
on a PC board as exciters. The structure is periodic, supported at discrete lines 
by transverse ribs. Figure 1 shows the electrostatic wall system along with the 
electrical hookup. Referring to figure 1, the output from the transformer T 
is biased at the center-tap and connected to the terminals A and B as indicated. 

The conducting membrane surface of width b and thickness h and isotropic 
properties (Young’s modulus E, density P, Poisson's ratio v) is supported 
structurally at separations of length S., where A « b; the membrane surface 
is electrically grounded. The electrodes are equally spaced from one another 
and at a separation H from the membrane surface, and they are connected 
alternately to terminals A or B. Each periodic bay of the conducting surface 
has a sealed cavity of volume V(H x £ x b) underneath it. 

The electrostatic wall model is designed to operate at a frequency range of 
200 Hz to 10 kHz. The largest dimension of the model is about 0.4 m, speci- 
fically for testing at the 7 in. x n in. tunnel facility at Langley. The 
largest nondimens ional speed is 
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wL 0.4 X 2ir X 10000 


« 1 


( 2 . 1 ) 


3 X 10^ 

where (o = 2Trf an<? c is the speed of light* Hence / the electric field 
max o 

between the electrodes and the membrane is qnasi-static- 

When the electrodes carry voltages ^ an electric field is set up between 
the membrane and electrodes and the membrane is subjected to a force field. 

The membrane deflects under this loading thereby altering electric field 
distribution* For the case where the electric field is only a weak function 
of the surface undulation, it is possible to uncouple the electric field into 
a primary field (field with no structural motion) and a perturbation secondary 
field which is dependent on the amplitude of the structinral motion. The deflec“ 
tion of the surface can also be split into a primairy deflection (under the loading 

due to the primary field) and a perturbation field, i.e. , 

E = Eq t El Cw,Eg) 

W = w{Eo) + wi(w{Eo)f E) (2-2) 

where Ei « Eq and hence wj « a)(Eo) 

3. STRUCTURAL ANALYSIS OP THE MEIffiRANE ON PERIODIC SUPPORTS 


3*1* Basic Approximations to structural Analysis 

We shall make the following assumptions with regard to the structiire under 
consideration: 

(i) The structure is a thin elastic membrane with isotropic properties 
(Young’s modulus E, density p and Poisson’s ratio = v) - 

(ii) The structure is rectangular, flat, and simply supported periodically 
at distances Z. Each periodic bay is identical with regard to the loading 
on it, etc- 

(iii) Beneath each bay of the structure is a cavity which is filled with 
an incompressible fluid* Hence any transverse motion of the structure should 
be so as not to decrease the overall volume of this cavity. 
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3,2* Juatificaiiion of the Basic Approximations 

We shall now examine the above assumptions for the case of a typical 
electrostatic wall as designed: 

(a) The elastomer used was a membrane of thickness between 2*5 ym to 25 ym- 
(Four thicknesses were actually used: 25 ym, 12,5 ym, 6,25 ym, and 2,5 ym) * 

The membrane was uniform and was an elastic material with v « 0 , 3 , p = 148 N/m^ , 
and E = 4 ?< 10® N/m^. In order to make the membrane conducting, it was 
aluminized on one side. The process of aluminizing did not change its uniformity 
nor its elastic properties in any considerable way, hence the validity of the 
assumption (i). 

(b) The membrane was usually stretched smooth and placed flat over supports 
which were nylon threads* The supports were spaced equally from one another. 

The nylon threads were cylindrical in shape and were glued to the membrane with 
a unifoina coating of epoj^y resin of negligible thickness. Iii all cases of 
construction the flatness of the surface was checked using an optical setup and 
was found to be extremely good. The excitation field on the membrane was 
obtained using a symmetrically arranged array of electrodes which were etched 
on a PC board. The electrodes had a flat geometry, were of identical thick- 
ness and breadth, and were equally spaced- The . alternate electrodes were 
connected to terminals A and E of a voltage source. Negligible current was 
drawn by the resistance of the wires, and hence a constant potential difference 
existed between terminals A and B. The electrostatic loading on the membrane 
was therefore identical between bays for the form of applied voltages and 

V (see section 4 for the analysis of the electrostatic loading) . Other 

B- , , 

forms of loading that might occur on the bays include fluid loadings due to 
static pressure differential appearing across the membrane (the static pressure 
difference can be held constant, or nulled using control valves for adjusting 
the back pressure of the cavities) and dynamic loadings induced by turbulence 
(if tested in a wind tunnel with flow over the membrane) and loadings due to 
fluid structure interaction, . The main body of the present analysis and 
esperiraental verification was for bench models where the fluid loadings were 
zero, and hence the fluid loading effects are not considered here. However, in 
a later section we include these cases and indicate how such cases can be 
included in analysis for an accinrate prediction of the ensuing structural motion. 
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(c) The cavity underneath the membrane was of uniform depth of about 
6 ym or 10 ym. For a membrane w’ith a width of about 0.2 m and periodic 
length of 0.01 m or less, the volume of the cavity is 2 10*^® m^ for the 

worst case. For a volume change of 0*01 percent, of the total volume by any 
form of motion of the structure/ the bay would have to suffer a loading 
of 1 kg/m^- Since the electrostatic forces acting on the membrane in each bay 
are orders of magnitudes lower in strength, such a change in volime is not 
possible during the motion of the membrane. In other words, the spring stiff- 
ness of the cavity for motions which tend to change the volume is many orders 
of magnitude greater than the stiffness (of the structure -h cavity) for 
motions where no such compression occiirs(i*e. the stiffness of the structure 
in a deep cavity or in vacuo). For the example considered above, the stiffness 
introduced by the cavity for motions where change in volume occurs is 


ki = AF/AV = 1 q 9 kg/m^. 


which is very large compared to the in vacuo stiffness of the stritture. 

Based on the above discussion it becomes apparent that assumptions (i ) , (ii) , 
and (iii) in section 2.1 are justifiable approximations to a structural 
analysis. 


3.3. Structiaral Response Theory for Transverse Motions 
of the Periodic Membrane 

The undamped structural response in transverse motion of a simply supported 
rectangular periodic membrane is given by 


, 3^ , 

ph + 

at2 


L12(1“V^) J\3x^ 3y^/ L 

Eh^ 


3^w 3^ 

+ N — - + N 


where the quantity 


X - 2 Y 9y2 xy 9x3y 


3x^ 


3^w~ [ 

3x3y 


(3.1) 


12(1 - 


is the flexural rigidity of the membrane, p^ is 


the total loading {in vacuo external load and the fluid loading) , W , N , DT 

X y xy 

are the midplane forces, and w is. the transverse motion of the membrane - 

To evaluate the midplane forces, one must have an idea of the nature of the 
loading on the structure. Based on a classification of the external loading we 
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can classify three distinct regions of structural response. To fix the above 
motion, let us use the folloisring nondimensionalizatiorx. 


Z 



Bj. - , X 
Eh S. 



Equation (3.1) in nondimens ional form is given as, 


(3.2) 


^4 

St2 L 


12(1 " \)^)J 


r r- 3^ - 

LaS^ a?^ J L ^ aP ^ ay^ 


+ N 


1 

z — ^ 


3y^ SxayJ 


- (3.3) 

For the case where the initial stretching forces applied at the edges of 

the membrane (i,e, the forces N W , N are all positive, the effect 

xo yo xyo 

of these tensile forces is to limit the amplitude of motion and thus raise the 
effective stiffness to transverse motion of the structure- When the initial 
fields ^yo' ^ 3^0 negative, the fields are compressive and instability 

of the structure may occur (i-e- buckling instability). 

When the initial compressive stresses are large enough to cause buckling, 
the rectangular membrane will no longer be flat, but will have a deformed shape 
under these compressive loads* With the presence of the cavity the initial 
shape of the membrane is very complicated possibly with local buckling within 
the bay at some points, etc* This situation is rather unpredictable because 
of precise knowledge of the initial tension or compression, especially for 
such flimsy structures* Since it is very difficult to apply compressive 
stresses while mounting, usually care is taken to keep the membrane with 
zero initial tension* The subsequent analysis and prediction are simplified 
by the assumption of zero initial tension* 

The following classification of structural response requirements is made 

with the assumption that the periodic bay is short compared to the width and 

the initial tension field is zero, i,e., h » Z and N = N = ET = 0* 

xo yo xyo 
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Regime I: small amplitude theory. - The region of applicability of this 

theory is for ~ = a/y <0.10? thus, this analysis is valid for load parameter 
n 

p < :_JE , response can be conducted using N = N = 

® 12 (1 - \)2) y3 ^ y 

- 0 in equa1n.on (3,3), The transverse motion is therefore decoupled from 
inplane motions for this case. 


Regime II: (intermediate) large response theory, - For this regime the 

midplane forces are no longer negligible* For analysis of this case. 

Von Karman*s theory is employed. For the type of flimsy materials we use^ 
even at very large amplitudes the structure is within elastic limits. Based 

on an order of magnitude analysis of terms such as u 2 compared to terms like 
2 

w we suggest that the intermediate large response theory we develop is 

A A- 

valid for ^ — or for load parameter p^ = 0,008. 


Regime XII: large amplitude theory (forp^ > 0,008- ) - For this region 

there are very few available methods of solution. The analysis should incor- 
porate large rotation effects. Equation (3.1) is not valid for such analysis. 
For static cases Reissner’s theory has been used with some regularity. The 
problem is a fully three-dimensional elasticity problem with all components 
of motion fully coupled with each other* Little progress has been made in this 
area over the years, Es^erimental studies in these amplitude ranges are few/ 
and most have been using plates which are fairly thick (compared to the thick- 
ness of 6 |im with which we are dealing) . Even before the small deformation 
implied in the theory for Region ii becomes invalidated, it happens more often 
than not that the tensile stresses under large amplitude motion go beyond 
tensile strength limits. The structure becomes locally plastic, and we no 
longer can use the isotropy assumptions implied in the analysis* 


For flimsy materials such as we use, the maximum tensile stresses induced 
in large anplitude static: motion are given by 



This stress should he belovi the tensile strength of the material we use, 
if assumptions of isotropy cire to be valid. For the periodic structure we are 
using in our analysis, the longitudinal intertial effects cire negligible as will 
be shown later in the analysis, and hence the large amplitude vibration problem can 
be viewed as a quasi-steady problem, especially while evaluating stresses. Hence 
the bound on the limits of an analysis such as that for Region II will definitely 
fail if 


t < '"ts 


where *^ts - tensile strength of the stmicture- 

Por the polyester film we use as the structural material the tensile strength 
^ts = 1.5 X 10^ N/m^ and E = 4 x 10® N/m^; the snallest wavelength used was 
X ~ 1.814 mm and the thinnest membrane had a thickness of h = 2.54 ym. From 
equation (3.4) we obtaine the maximum limit on w/h beyond which the structure 
will be at least locally plastic as 


w 

h 


4 


Ts*4*(l - v^ ) _ 2 


E X 


Try 




(1 -V^)'^Ts 
E 


84.0 


(3.5) 


We had given the limit for validity of analysis for Region H by consideration 
of the small deformation approximation implied in it as 


w 

h 


0.03 

Y 


21.5 



(3.6) 


for a = 1.4 X lo 
we are designing 
analysis for the 


a typical case. Hence, it is feasible using the structure 
to study the limits of validity of the moderate large amplitude 
first time in our knowledge. 
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3,4. Structural Responae Thsory for (Intermediate) 

Large Responses of the Periodic Structure 

Equation (3,3) in nondiitiensional form represents the governing equation for 
transverse motion of tt*e structure. We make assumptions (i) , (ii) and (iii) 
of section 3.1 in our analysis? i.e. b » £ and 

fj'\f dx dy = 0 

V 

furthermore, we assinae that the excitation field is of the form 


f p sin dx =5 0 for m = odd numbers (3*7) 

0 ^ 

where is the loading function in each. bay. The excitation field that will 

be derived in the next section can be shown to be of the above form for a given 
form of voltages at teanninals A and B of the electrostatic setup. 

For a simply supported rectangular bay, the deflection shape can therefore 
be only of the form 



ran 


w 

ran 


sin 


2m7TX 

£ 


. sin 


nTTT 

b 


(3.8) 


Furthermore, because of uniformity of the electrostatic field in the y direction 
the choice of the deflection shape is restricted to 

w =25Zw sin sin (2n - 1) ^ (3-9) 

. m n ran £ b 

or, in nondimensional form, 

a =EE a sin(2m7Tx) sdLn(2n - 1) (iry3) 
m. n ran 

where S = £/b (3.10) 

defining: 

m = 2iinr (3.11) 

(coht^d) 
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yields 


= (2n « Dgir 


a=EE 


(t) * “ - . - - 

a san m x sxn 0-y 
mn 


(3,11) 

(concl’d) 


3.4,1. Midplane forces. - During deflection under transverse load the 

midplane suffers in plane strains E , E E - i^ssuming that the external 

XX YY/ xy ^ 

load parameter extends through Region II, v;e define the midplane forces as 


N 


X 


N =: 

y 


(1 - v2) 


(E + 

(1 - v") V ^ 


V E 


yy> 


V E 


XX , 


N 

xy 


Eh 

2(1 + V) 



(3.12) 


In terms of the nondimensionalization adopted in equation (3.2) ; 


N 

X 


= E — + V E~ 

1 _ l_-xx yyj 

1 „ ^2 L yy xxj 


N 


B— 


xy 2 (1 + v) xy 


(3.13) 


The raidplane strains are related to the components of deformation as 



, ^ ^ ^ ^ ^ ^ 

^ By 8x 3x 3y 


(3.14) 
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where u = u/^, v = v/5- and Ur v are the dimensional in~plane displacements. 
Before setting up the dynamic equations of in-plane motion we introduce the 
following consideration: 

The frequency of excitation for the electrostatic wall is well below the 
longitudinal natural frequency of the structure. For this case the in-plane 
motions are stiffness dominated. Hence, the effect of longitudinal inertia can 
be neglected. 

The above requirement is satisfied if the reduced frequency « 1 

(but a more relaxed criterion can be < 0.10) , where f is the operating 

frequency, c^j^ is the shear wave speiS, and L is the wavelength of the 

primary motion. For the electrostatic model L = A, the maximum frequency of 

operation f =* 10 kHz, and c^^^ 1690 m/sec. For the largest wavelength 

model that we considered, A = 7 mm. Hence, = 0.04 < 0.1. Thus the 

^sh 

neglect of longitudinal inertia in all the siibsequent analysis is well justified. 
The dynamic equations of in-plane motion are now given as 



3x 8y 


3N- 3N— 

—L + — 32 := 0 

3y 3x 


(3.15) 


Substituting the relations given in equations (3.13) and (3,14) we get 


3^u 

3x2 


1 - V S^u . 1 + V S^v 


Sy-" 


3x3y 


+ F-(x,y,t) - 0 


3^v . 1 - V 3^v . 1 + V . 

' ' ■ “ ”T ■ -j- ” ■ ■ ■ ■ ■ ^ 

3y2 ^ 8x2 ^ 3x3y 


p-(x,y.t) 


= 0 


(3.16) 


where 
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9a Ta^a ^ 1 - V a^-a l ^ 1 -f v 3a 3^a 

3x Lsx^ ^ 3y^J ^ 3y 3x3y 

3a r 3^a ^ 1 - v 3^a ~j ^ 1 + v 9a 3^a 

3y Lay^ ^ 3x^J ^ 9x 3x3y (3.17) 

For the deflection, shape, we use equation (3.11). Substituting in ec[uatiQn (3.17) 
we get expressions for F- and P- as: 

= jEEEi: L<^ * 3>^] ®i’y] 

• [3 ■ 5 >*] 

• “=['% * 4 * 

- sin[(I + j)S] cos[(ejj - %)y][j{i^ + s| - 6j 6j)] 

- sin [(I - JlS] cos [(Be - %)5][5{P + ^ S| - Bj Bj)] 


E’-(x,y,t) = 


F-(X/yrt) = 


= jEEEE j sin [(6- + e-)y] cos [(I + j)x] 

■ ^ 2 ^ 2 ^ ■*■ J °°® 

* 2 2 ^ ■*■ E^^k “ E^^ 

* E^f ^ ~2 ^ ^ 2 ^ E^^k “ E^i - j)xj 


[Bj( 6| + i !)]) 


(3.19) 
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From the Icnown form of F- and F- one is now able to proceed to 
solve u, V. 


3,4.2. Edge conditions , ^ The edge conditions describe the state of 
fixity of the structure for in-plane motions along the edges* For the simply 
supported periodic structure the edge restraint is 

^(OrY/T) = u{l,y,x) 0 

v(x,o,T) = v(x,|-,x) = 0 (3.20) 

3.5. Solution of the In-Plane Motions 


Since equations (3.16) are time-dependent equations, a transient solution 
which depends on the initial condition will exist along with the solution of 
equations (3.16) . This particular solution will be representative of the initial 
displacement fields that may exist (e.g. , if initially the membrane is stretched 
with uniform tension T, the initial displacement field will be nontrivial) . 

For the case of zero initial tension, hhe particular solution takes the 
trivial form: 

u^(x,y,T) = v^(x,y,T) = 0 (3.21) 

For the solution of equation 3.16 we take the displacement fields as. 




sin 1 

(i + j)x 1 cos 1 (0- 



- j)x 

1 COS ^ 

((eg + 0j)y] 

1 


+ j)x 

1 cos 1 

(Bj - Sj)y 1 

1 


- j)x 

} cos 1 

1 <9k - 1 

Q 

(3.22) 
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V =EEEE coajd + j)x} sinjte- +e^)y} 

■*^ cosjd - jl^J Sin|(3- +0^)y} 




Substituting for the derivatives of u, v, P-, P- in equation (3.16) and 
equating like terms one obtains a set of alg^raic equations for the coefficients 
■^ijkA' 5 ive one set of these equations as: 

[(^ *i) J ^ [H^ j) (^S ^ 65)] 


ip ^ e| * sj Si 


(3.24) 


^ijk£ (^k ^ ^l)(^ ■*• ^)] ®ijkil[(pk H) ^ ^ 


4 ^ik 


jA[^£ 


, 1 — V » 1 + V T T 

6 j + 1 1 3 


(3.25) 


and similar equations connecting other coefficients. Finally the coefficients 
are given as: 
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For i / j / 


* = + 


a., a.- 
ik 


ijki, 2(1 - V) „ 1)2 + (P_ + g_)2j 


4 * 8 j] {(®s ^ (^ - 5 )' 


+ (^ - 5 )(®k * ®j) (4 + 


1 - V 1 + V T ~ 


i — ' ■ ' ' X. 

2 2 


J)] 


(3.27) 


when i = j, - 0. 


C..,.„ = - 


a * 1 S * fl 
XK 


ijkS, 2 (1 - \i) 


L(i+ j)2 + (3- - 0^)2J 

.|j[i" -H s| - ej 6j][(6j - Bj/ + ^i^(i j)'] 


+ ifi (l + 5 )(s5 - 6j) 8j (b| + 




1 + V T t\\ 

- 1 -^Vj 


(3.28) 


and whenever i = j, ^ other cases 


a,, a. - 

1 - xk 

ijkJl " 2(1 - V) 


n 2 


(I - j)^ + (3j^ - e^)^J 


^ 1 - V 2 1 ±^_ p p 


i) l^k 


k •*• 2 


-V — 

-2.-3 J 


•t(r 

. (l - j) 6- (e- - e-)(4 + - i-p I j)] ,3.2„ 
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— Vf 


L(i + j)^ + (aj^ + a^)^j 


(®k + sj) * 5 ) j { e| + ej 6j) 


a., a. r 

P _ ajc 2% 1 ; 

iJWL “ " 2(1 -V) [,J . J,2 + (6j+ Sj)2 


•{H^ (sk * Sj) (j - 5 r 


j) j I i + ^ 2 


1 + V 1 

2 ^ic ^Jlj 


[i* - i) 


" * (sj + %i 


0c + 


X - \> 1 + \j T rT 

1 - 1 — ^ 3 J 3 ^ 


(3.31) 


= 0 for k = J, and 


ijka 2(1 - 


ik^j£ 1 

^ L(i + j)2 + (3j^ - &-)K 


(^k - Bj) (^ * 5 ) 3 [i^ * S^ 0 | - 

Bj[(l + if + - 0j)2] [0I + - 


1 + \> 

~2 8k ej 


. 1 + V T T 

1 3 


(3.32) 


and H. ,, „ = 0 if k = £ or else 
i]k£ 


a. a . 
ik i£ 


xjk£ 2(1 - V) [(j _ jj2 ^ _ g_j2 


■[H^ k - 82 ) (J - 5 ) 5 ^ i 8 J J 

- 82 [^ - i)" k - 82)8 H4][si 


1 - V 1 + V T “ ( 

(3.33) 


L6 


Having obtained the coefficients one is now able to evaluate the laidplane 
forces through the help of . equations (3,13) and (3,14) . Thus 


N = 

X 


(1 - V2) 


h T 


^ + V E.j 


ijka 


ijka 


(^k 


(l j + V g- j cos(i + j^x cos^Bj^ •{• Bj^y 

(®ijU ■ i) + '■ '’iJM (®S ■" ®s) ■ (i 5 - V Sj 6j)) 


COS ' 


(i - j); =os(6jj + 

\jM (»k - ®x) i ^ 8l)l 


cos 


(I + j)x cos(bj^ - Bj)y + jcos(i - j)s cos(b|^ - Bj)y) 


{ °ijk£ " ^) ■J' ®ijkA (®ic “ ^l) ' 8 ^ 


N = 


j'’ \jki + i) + (®k * ®») 


Cl - 


- (Sfc ^ j) 1 “= + j)* =“ ki ^ 

*{'■ - i) + ^ + fi!^ . J 5 „jj 

• C=,s(i -j)icOs(^ +S5)j + |vCyj^ (i + j) 


(3.34) 


V I j - B- 6 J I cos (i + j)x cos (b- - B^)y 


(3.35) 
(cent *d) 
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J 


“UM (i - =ij« (Ss - .^l) * ['■ i J \ 


cos (JL 


- j)5 cos(gj^ - e^)y 


N — — 

xs' 2(1 + V) 


'^[{-^ - k*^ 


- =iji« (^ + i) ) + J)^ =“(^ + 


B. „ 

X]k£ 


(»E *6^ ■" "'ijk£ (i - j) ■" Be j } =in (l - j) 


“(Sk * Bj)5 +(- Cykj k ■ (^ * ^) 


xk 3£ 


ej j} sin(i + j); sin(Bj - B^f - (bj - 6^ 


* “ijkl (^ - i) ♦ be 5 } - j): sli.(Bg - 


Having obtained the terms N , N , N one can now solve ecniation (3.1) . 
^ x' y xy 

3.6. Solution of the Transverse Motion 

The governing equation (3.3) can be solved as modal equations in time. 
The form of a is given by equation (3.11) as; 

a =T’y' a ^ sin m x sin (6a y) 
xnxi 


The Raleigh-Ritz approximation to equation (3.3) is obtained by multiplying 
equation (3.3) by sin p x sin g~ y and integrating over the x - y space. 
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Defining: 


\jkA 


ijU 


ijkjl 


and 


3 < 2 > = 

ijkl 


(i + i) 

^ijka \jkii 


- -^ (IJ ^ V aj sj)] 

(3.37) 


(i - j) + V 

(es 


(3.38) 

[*^iDkJ!. 

% 

(l + j) + V Gy,^^ 

(«k ■ 8^ 

+ (i 3 - V aj 8^] 

(3.39) 

[°ijkP, 

(i - j) + V Hy,^^ 

(«£ - 

|.^(iJ-vaj8j)] 

(3.40) 

^ijk£ 


(pE* ej) 

- (a, 8j . V U) 

(3.41) 

®ijk£ 

(" - 5) 

(»E - 


(3.42) 


^ 5) + Gyj^j 

(®E ■ ^i) 

+ — gj- (v 1 : - Bj ej) 

(3.43) 

°ijkil 


(«E - «e) 


(3.44) 


and 


=«L = 5 5 - \jicj (be Sj) - (l 4- j) 

' - [«ijM (sE 6 e) '■«!« (f - j) Bj- j] 


ijk 


(3.45) 


(3.46) 
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(3.47) 



and also the following- parameters 
Rl = i -r. j + p 
R2 = i + j + r 
R3 = i + j - p 
R4 = i t p - j 
R5 = i + r - j 
Rg = i + p + r 
R7 = j + p + r 

and 

Sx = k + S, + g 

52 = k t A + s 

53 = k+ s+ q 

54 = £ + s + g 

55 = k + “ g 


(3.48) 


(3.49) 


(3.50) 

(coat'd) 
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1 


Sg — Jt — A + s 


Sy = k - s + g 

and the Kronecker delta fiinction 
^rs = 0 r s 


(3.50) 

(concl'd) 


rs = 1 r = 1 


(3.51) 


the Rayleigh-Rits^ approximation equation (3.3) is ■ Hipn 


d'^a 


_EH. 


+ K a + 




^^2 P'2 pg i j k £ r s 4(1 - v^) 




* fil 


?1 + §2 P£ + 83 P 3 + 64 




, pq 
‘Gen 


(3.52) 


where 


K = ^ ! — 

12(1 - v^) 


(? . 


(3.53) 


Si = 


a a 
ij k£ 


p^ 

Ijkil “ q ijkJl 


P2 = ^ a«2) 1 

^ij^£ L ®q ^ijkij 


P3 = 


^ij^k£ ^ 


4 ftt3) 

^ ^ijfci ^ ^q ^ijk£ 


Plf = 


^ij^£ 


"u2a^^^ + =(4) 

^ ijki ®q ^ijk£ 


(3.54) 


and 
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5t = (e +6-6 \ (s +6-6 ^ 

V Rlr R2P Rsry \ Sis S2q Sss/ 

52 = + 6 ^ - 5 ^ ') / 5 „ + 5 ^ - 6 ^ . - 6 ^ 

\ Sis Saq S 5 S/ \ Ri).r Rsp Rg^ B.7V 

53 = “ 5„ ) + 5^ - 5^ - 5^ 

\Rir R 2 P Rar/ V SyA Sgq 33/1 Sijk/ 

54 = + 6 - 6 . . - 6 A /6 +6 - 6 - $ ,] 

^ V R4r Rgp R6 j R71/ \ S7il Seq S3JI Sifk/ 


and 



(x,y) sin (p x} sin {g- y) d- 

q X 



(3.55) 


(3.56) 


!The set of ordinary differential equations is highly nonlinear because of the 
presence of terms such as ^ij^A ^pq. The coupled set of such equations can 
be solved using numerical techniques. Itoalytical solutions might be impossible 
except in the simplest case/ e-g. a one mode solution such as i = j = k — A = 
p = q = 1. When the electrostatic wall model is subjected to fluidic forces such 
as will occur if used for drag reduction studies, a complete analysis of the 
problem as formulated above with a few normal modes in the deflection shape might 
be necessary. This is especially so in flutter regimes of structural motion. In 
such cases the generalized loading will consist of three main contributions: 

(a) random turbulent pressure loading which can be insignificant unless flow speeds 
are fairly .large and the membrane thickness very small; (b) wall pressure loading 
due to the interaction of a pulsating boundary in a turbulent boundary layer flow; 
the magnitude of this loading can be significant under a variety of circumstances 
including the flutter mode instability case; and (c) the electrostatic wall loading 
due to the primary excitation through electrostatic attraction, forces . 

In an actual experiment with flow over the structure such as in a windtunnel 
other loadings ■ might also be present. The most common of these are pressure 
gradients in tne tunnel and static pressure differentials across the membrane. An 
accurate prediction of . the structural response depends largely on identifying all 
these influences,. 
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4. AiaLYSIS OF THS ELECTBXC FIELD 

4.1. Determination of the Primary Electric Field 

The schejnatic of the electrostatic system is showi in figure la. The output 
from the transformer T is biased at the center tap and connected to terminals A 
and B as indicated. An array of electrodes which is etched on a printed circuit 
board and coated by degassed epo:^ cement or lacquer forms one part of the 
electrostatic wall system. A grounded conducting sheet (membrane) which is struc'-* 
turally supported periodically at distances apart and separated from the 
electrodes at a height “H" forms with these electrodes a capacitive network. 

The membrane is free to move transversely; when alternate electrodes are connected 
to the terminals A and B respectively the membrane is subjected to an electro- 
static force distribution and deflects under these forces. The space between 
the membrane and the electrodes is filled with air and sealed* To obtain maximum 
force for given terminal voltages, the distance between the membrane and the 
electrodes should be as small as possible. However, as the distance between the 
membrane and the electrodes becomes smaller and smaller the applied voltages in 
the terminals must be reduced lest a breakdown of the electrostatic wall system 
due to arcing will occur. The breakdown potential for a given separation in air 
is given by Paschen^s law- For air at NTP the breakdown rms potential gradient 

is 3.1 kV/itun. If instead of air a fluid such as SFe is used, the breakdown voltage 
gradient can be raised by a factor of two. 

In figure 2 we show a typical bay of the electrostatic wall configuration. 

The electric field between the electrodes and the membrane can be obtained by 
solving in the domain 

v^v =0 (4.1) 

with the boundary conditions (as indicated in fig. 2) . 

V = 0 

V =Vo + Vi(x,H) sin nt 

V = Vo 2/H (4-2) 


z - 0 

z = E 

X = 0 , 2 , 


In figure 3, Vj (x,H) and its Fourier series representation is shown. 
From the figure it can be seen that 
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(4.3) 


Vi(x,H) = 1.1463 Vj sin 2iTy sin flfc 


is a good approximation. 


The solution of equation (4-1) with the boundary conditions given in equation 
(4.2) is rather straightforward and can be written as 


V Z „ sinh (2Trf- ) 

V = — + 1.1463 Vi sin sin fit ^ 

“ ^ sinh 2 tt f 

Z 


(4.4) 


The electric field in the dielectric between the electrodes and the membrane 


S „is thus 


E = E = o 
X y 

V Vi cosh[l2z] 

E^ = - ^ 7.2024 - sin — sin fit — 7 ^. 

smh( — ) 


(4.5) 


Force density due to primary field. The membrane is stabjected to a force 


p® - — ^ I 

o 2 'z ^ 0 


(4.6) 


Thus the generalized electrostatic force is 


= 0 P ^ ^ 

Gen 


and for P = 1 


ppq L- ro>935 £ V n 

“■> “ (20 - 1 ) .ii.h(^) L -I 


sm 0 ) T 


(4.7) 


Optimal values for Vq, There is a limit to the maximum voltages 

Vq/ permissible for a given conf igviration , which is determined by the imis 

voltage gradient to be less than the breakdown value (3.1 kV/mm for air at NTP) . 
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From equation (4^5) the maxiimim ntis electric field value is given as 



£ 3.1 Mv/m (4.8) 

In order to obtain the optimal values of Vq, for maximum generalized 

forces^ one must minimise the generalized forces subject to the constraint 
equation (4.8) . 

4-2. Requirements, Ratings > etc, of the Electrical 
Networks for Electrostatic Wall Configuration 

In figure 4a is shown the schematic of the electrical setup. Power is 
drawn from a signal generator stepped through a power amplifier and fed into a 
transformer* The transformer output is fed into the exciter for the electrostatic 
configuration. Basically the electrostatic wall comprising the membrane and the 
terminals act as a capacitive load- The capacitance can be measured using a 
capacity meter available commercially- It is also fairly straightforward to 
theoretically model this electrostatic wall as a discrete capacitive network and 
evaluate the capacitances of the wall as done in reference 6, 

For satisfactory perfoinnance of the electrostatic wall system the fallowing 
points should be kept in mind: 

(i) The electrostatic wall should be operated below the breakdown voltage 
levels to avoid arcing and sparking and consequent degradation of the terminals 
and burnout or charring of the membrane surface* 

(ii) Collapse of the membrane into the cavity and consequent shorting of 
the electrical system should at all costs be avoided lest the transformer or 
power amplifier be damaged- 

(iii) Any L-C oscillation due to the loading of the system should be avoided. 
This is accomplished by designing the outages of the units such that the external 
load is well within the operating load for the system* 
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4*3. Evaluation of the Maximum Peimiissible Capacitive 
Load for the Electrical System 

The ratings of the transformer are 

Turns ratio a = N 2 /N 1 = 40 

Resistance on the primaary side - 0-5 Q 

inductive reactance on the priiuairy side Xi = 8 @ 1 kHz 

Resistance on the secondary side R 2 = 2 

Inductive reactance on the secondary side X 2 — 50 § 1' kHz 

External load is capacitive with a reactance X @1 kHz “ 

^ C ^(liC 

In figure 4b the equivalent circuit is shown- It is assumed that the 
magnetizing current is zero, an assumption which is quite valid (usually it is 
5 to 10 percent of the primary current) * With this assumption the reactance 


From figure 4: 

The equivalent impedance on the 
primary side = Zeq 


= j X + Ri + j — : 
Ll ^2 


L2 R2 
,2 =2 






The primaary current 


V/z = - /-0 
eq z ^ 


The load current 


l 2 = Ii/a = — ^-6 


The circuit becomes a resonant circuit when the reactance is zero. 


i-e- X^ t X, = 

^2 t.ca2 


The operating range of the electrostatic wall for experiments lies between 
200 Hz to 2 kHz. In order to protect the power amplifier the resonant frequency 
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which can be determined from equation (4.12) should be well above this range. 
Fixing the resonant frequency f to be around 3000 Hz we obtain the condition 
that capacitance of the electrostatic wall be 


c < 


2TTf 


■Xjf 

1000 


a^ + 


X2f]' 

1000_ 


280 pf 


With a proper choice of equipment (i.e. low reactance transformers and 
power amplifiers) , the size and capacity of the electrostatic wall can be 
enhanced considerably for the same range of operating frequencies. The fore- 
going discussion makes it imperative that when different models require 
laborato3:y testing the capacitance of the model be measured to see whether 
there is any limitation on operating frequency imposed by the choice and 
availability of power equipments. 


5. THEOEETICAL ANAIiYSIS AND EXPERIMENTAL VERIFICATION OF THE 
PERFORMAJ.^CE OF THE ELECTROSTATIC WALL 


5.1. Analysis of Structural Motion 

We consider the case where the generalized force is due primarily to the 
electrostatic field set up in the model. The structraral motion resulting can 
be studied by solving equation (3.52) with the expression for p^^ provided 
by equation (4.7). The deflection shape a is given in equation (3.11) as 

a a sin m x sin S- y 

m n mn n 

We consider the case where 6 = i/b « 1, corresponding to the situation of the 
width of the bays b = 50 fi. which was true for most of the wall to be built for 
tunnel testing. 

From the nature of the generalized load given in eqxiation (4.7) / 

Gen 

i.e. 


P 


pq _ 

Gen 


0 P > 1 
P(q) p = 1 
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it is obvious that the deflection shape can be approximated as 


a « a sin 2 Trx sin 6 - y 
n m n 


(5,1) 


It is not at once clear how far the summation implied in equation (5.1) must 

be carried out. However^ the natiare of the generalised forces (being 

1 Gan 

proportional to -* 1 ) ^ suggests that at best the first two modes are important. 
The equation system for two mode series is of the form 


“ + an [kii] + Gh + g \^ 


+ ^22 [ k 22 ] + G ^2 ^2 + ®21 ^11 ^12 + ^23 ^12 = 


(5.2) 


In table 1 the coefficients ch/ Gi 2 r G 13 , etc. are all tabulated 

for arbitrary 3, using a program called MACSYMA available through AKpa., NBS 
or MIT network. 


S . 1 . 1 . Initial approximations . - We shall examine the case of 3 « 1 
and assume that the percentage of the second normal mode in the solution to the 
problem is entirely insignificant in order to obtain an idea of the sensitivity 
of the solution to tanmcation of the series. Setting 3 = 0 and taking only the 
first of the equation set in eqiaation ( 5 . 2 ) one obtains 




+ kii an t Gh a-?i = 


(5.3) 


From equation (4.7) ^Gen. ~ 

^ 0.935 e Vq V 1 I 

where ^10 2kK EyH2 


(5.4a) 


(5.4b) 


Then terms kn and in equation (5.3) have (for 3 = 0) values 


2S 


^11 = T ' ‘^ll “ ~ ^ 

12(1 - V*) (1 “ V^) 


(5.4c) 


Equation (5.3) can be rewritten using the fallowing transformation: 


B = I - n Tr u = a.ii .^1 , = 


Q 


kii 


and P2 = ^XO- 


/^U 


(5.5) 


as 


n 

— ~ + $1 u + u^ = P2 cos S 
d6^ 


(5.6) 


5.1.2. Harmonic solution of equation (5-6). - Seeking harmonic solutions 
for U/ equation (5.6) leads to a simplified cubic ec[uation. 


(01-1) u + ^ u^ = P 2 


(5.7) 


from which u, and hence a, can be obtained. 

5.2. Theoretical and Experimental Results 

Theoretical results for two configurations of the electrostatic wall are 
presented in figures 5 through 8. The membrane material used (commercially 
named Mylar) has the following properties: 

E = 3.5 X 10® kgf/m^ 


p = 138 


kqf sec^ 


m^ 


V = 0.3 


(5.8) 
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The air gap between the membranes in two dixferent configurations were 
H - 0,305 mm and H = 0,127 mm. Since the electrostatic forces are inversely 
proportional to gap height H, the smaller gap provides a larger force on 
the structure for a given voltage. However, the breakdown potential for the 
system is smaller when the gap is smaller. Hence the electrostatic forces 
available in both these configurations at the optimal voltage are of the same 
order. However the smaller air gap has the advantage that one can work with 
lower voltage levels. 

In figures 5 and 6 the theoretical values are compared to the bench test 
results- It is surprising to observe such an excellent agreement between the 
test values and theoretical prediction. The experimental measurements of 
dynamic surface motion were obtained using an optical setup. Reference 10 
describes in detail the measurement technique. The a.c. voltage applied to 
the electrodes was at 300 Hz for these experiments. 

Figures 7 and 8 show the frequency response of the structure at various 
excitation levels. The nature of the backbone curve is indicated by dotted 
lines in these figures. At low excitation levels where w/h « 1, there is a 
steep increase in amplitude levels at resonance. At higher levels of excitation 
resonance has little effect on amplitude level. In the theoretical analysis 
structural damping was taken to be negligible* Since damping primarily affects 
the near resonance amplitude levels, it is not necessary to include structural 
damping in the response studies for . the present case as can be seen from the 
nature of the response curves. The reason for this is very clear; viz, the 
test structures have resonance frequencies well above the excitation frequency 
for load values of interest (large amplitude case). Furthermore at these load 
values the nonlinear stiffnesses themselves act as a damper or delimiter on 
the amplitude levels. When the structural nonlinearity is of the soft spring type 
it might well be important to include the dampdjig in the analysis since natural 
frequencies for this case are being constantly shifted dovmward from the low 
excitation case (w/h « 1) with larger exccitation levels. 

5.3. Examination of the Perturbation Field 
Due to Primary Motion 

When the membrane is set into motion by the action of an electrostatic 
force field, the electrical field configuration changes and consequently the 
excitation field itself must be reexamined. 
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To analyze the effect, of the surface motion on the electric field configu- 
ration an orthogonal curvilinear coordinate system can be used. Since, as 
far as the electric field is concerned, the structural motions have a station- 
ary behavior, a quasi-rigid coordinate mapping is required. Curvilinear 
coordinates such as the one to be described have been used by Benjamin (ref. 7) 
and Chang (ref. 8 ) for flow problems. The first order curvilinear orthogonal 
map is shown in figure 9 . The amplitude of the wavy surface is taken to 
be ( 0.1 h); Y]_ = 0.0 curve in the figure represents the approximation to the 

actual surface. Since — « 1, the variation in the chordwise directions have 

b 

been neglected just as in section 4. The mappings for this configuration are 


XX = X - a e~^^3 CQ 3 sin Qt 

yx = z - a sin 2ttx sin J2t (5.9) 

X - y coordinates can be written as 

X = XX + a e cos 2‘irxx sin J2t 


z - yx i 3- 6 2irxx sin S3t 


(5.10) 


The solution of equation (4.1) with the boundaiy conditions should be attempted. 
We write the governing equation and the boundary conditions in the mapped system 
as 


a2v 

— — + = 0 


3xx^ 

9yi^ 


Yl = 

0 

V = 0 

yi = 

(HA) 

0 

> 

(1 

> 

XX = 

0, 1 

■ — 1 

> 

il 

> 


(5.11) 


(5.12) 
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The solution of the above problem yields 


^2 H 


^ \T' 


sin 2ttx sin Qt 


+ 2va (cos 27TX __ qqq 4ttx) sin 


V - 0 
X 


The correction to tlie electric field distribution is therefore of the order of 
(w/ii) times its original value. For the electrostatic walls we are testing 
the magnitude of this correction is less than five percent when the amplitude 
of the surface motion is 0.025 mm at a spacing of the bay of 1.S14 mm. In all 

the tests that we carried out the perturbation field was negligible since the 
excitation field could barely drive the system to these values of amplitude. 

An error of 15 percent in the evaluation of the electric field would usually 
generate an error of 10 percent in the amplitude prediction, which must be 
borne in mind when accuracies of that order are required* 

When the deflections of the surface are quite large the perturbed field 
can be evaluated using a more refined (second order) coordinate system such as 
used by Chang {ref. 8). We show the mapped system in figure lo and give the 
transformations below as 


XI = X - a cos 2irx sin Sit - ira^ sin 4irx sin Sit 


YX = z ^ a e sin Zirx sin Sit + 7ra^ 


cos 4ttx 


^ 2 C 1 4- a cos 2 ttxi “ eT^'^l s±n 4'irxi 


Yl -f a e 2ttxi ira^ 


47TXX “f 1^ 


(5.15) 
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5.4* Inclusion of Aerodynamic Forces for 
. the Electrostatic Wall 

The electrostatic walls were purposefully designed as small wavelength 

configurations. Hence, static divergence will not be achieved in any of the 

electrostatic walls under operating condition {15 m/s - 50 m/s) in the 7 in. x 

11 in. tunnel. In reference 1, an inviscid analysis has been carried out for 

aero 

the evaluation of the aerodynamically induced pressTore p [eq. (39) , p. 16 of 
ref- 1] due to structural motion. 

We give this expression as 

q2 

aero „ . “ . - — 

p - 2?rp<arr> an sin 2tfx sin Tiy sin ojt. (5.16) 

Vl - 

Figure 11 shows the effect of including aerodynamic loads in the analysis 
for the case of flow over the membrane. Our contention is fairly clear; there is 
no need for any sophisticated analysis of aerodynamic effects for these configu- 
rations since even the inclusion of the inviscid values which overestimate 
these magnitudes increases the levels of amplitude very slightly and that increase 
is only over a very narrow window. The’a.c. voltage applied to the electrodes 
in this example is at 300 Hz. If the* a.c. voltage was applied at very much higher 
frequencies the effect of the aerodynamic load- on the motion will be enhanced 
since the effective nonlinear stiffness would be considerably lowered. 

We again point out that the analysis which took into consideration the 
aerodynamic load was carried out under the prior known fact that the wavelength 
of the structures was much too short to cause static divergence; hence tlie flow- 
structure interaction problem is not of an eigenvalue type but merely a forced 
response problem. 

6 . CONCLUDING REMARKS 

A unified theory for an electrostatically driven active wall system has 
been presented. The electrostatic wall system, is capable of producing deflections 
of many orders of thicknesses. Consequently a large intermediate response theory 
has been used for analysis. The theoretical analyses; are compared with bench test 
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results which show excellent agreements between the two* The case of an elec-* 
trostatic wall vibrating in the windtunnel is considered, and the theoretical 
predictions under simulated flow conditions indicate that the aerodynamic 
effects are negligible* The perturbation effects to the electric field due to 
the structural motions are also considered and shown to be negligible for the 
test es^eriments * 
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Table 1. Stiffness coefficients. 
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Figure 6. Theoretical response curve and test results 
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Figure 7. Frequency response curve 
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